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ABSTRACT 


Numerical solutions for radiative transport in a class of anisotropically scat- 
tering materials are presented. Conditions for convergence and divergence of the 
iterative method are given and supported by computed results. The relation of two 
flux theories to the equation of radiative transfer for Isotropic scattering is discussed. 
The adequacy of the two flux approach for the reflectance, radiative flux and 
radiative flux divergence of highly scattering media is evaluated with respect to 
solutions of the radiative transfer equation. 
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INTRODUCTION 

Radiative transport theories involving multiple scattering play an important role 
in the engineering analysis and simulation of the performance of diafhermanous materials. 
Highly scattering dielectric materials have been proposed, for instance, for the entry 
heat protection of planetary probes [1]. Other applications include the evaluation of 
the reflectance of condensed deposits on cryogenic storage tanks [2] and the engineering 
analysis of opacity in the paint and paper industry [3]. Regardless of the application, 
engineering approximations of radiative transport continue to be employed despite the 
increasing ease of computer modeling. Accordingly, this work attempts to provide a 
unified discussion of the radiative characteristics of anisotropically scattering volume 
reflectors through new numerical solutions of the equations of radiative transfer, to 
evaluate several approximations for parameters of importance in thermal analysis, and 
to clarify the ties which exist between the equation of radiative transport and the 
approximate theories. 

. Of primary interest here ?s the investigation of the radiative energy transport 
characteristics of volume -reflecting heat shields for hypervelocity entry. In this 
application a highly scattering, weakly absorbing dielectric material mounted on an 
opaque, structural substrate is irradiated by an intense, diffuse radiative flux. Low 
absorption allows the incident radiation to penetrate in depth with limited conversion 
of radiative to thermal energy. In addition, high scattering provides an internal 
mechanism for spatially distributed reversal of the direction of the incident radiation. 

In this way, dielectrics are capable of reflecting intense incident radiative fluxes 
efficiently while convective heating is also handled efficiently through ablation mecha- 


nisms. 
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Multiple scattering studies including effects of anisotropy have focused in the 
past on the reflectance, transmittance, and emittance of plane "parallel layers exposed 
to parallel incident irradiation [4,5,6]. Solutions of the transfer equation presented 
here assume a diffuse, incident intensity field appropriate to the entry vehicle appli" 
cation. With this orientation, emphasis is placed on thermal analysis parameters and 
their approximation. 

The approximations considered for isotropic scattering stem from the Schuster- 
Schwarzschild [7] and Kubelka-Munk [8] two-flux models. Exact solutions of the two- 
flux equations as well as approximations [9] and a reinterpretation [10] are evaluated 
in relation to radiative transfer equation solutions. An earlier evaluation of the 
Schuster-Schwarzschild model [1 1] considered coupled radiative and conductive energy 
transport in a broader context whereas the present study focuses on the uncoupled, high 
albedo problem. 

ANALYSIS 

Transfer Equation Formulation 

The Equation of Radiative Transfer is considered to govern the transport of 
radiation in dielectric media. Polarization effects are known to be negligible In 
thermal problems involving multiple scattering [12], and are therefore not considered. 
A diffuse, gray radiative flux is assumed. The medium is assumed to be isothermal, 
to be gray, to have unity index of refraction and to have uniform absorption and 
scattering coefficients. The rear boundary is opaque and reflects specularly with 
constant reflectance Rg. 

In order to evaluate the influence of anisotropic scattering, a phase function 
of the form P ( © ) = w ( 1 + x cos 0) is assumed [6] . Here to is the albedo, 0 is the 
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angle between incoming and scattered beams at a point in the medium; and x is an 

anisotropy parameter which provides backward scattering (-l.O^x *^0), isotropic 

scattering (x = 0), and forward scattering (0 < x ^ 1.0). 

-The intensity field is azimuthally symmetric in this plane-para! lei case and 

may be represented by 1 { T , M’) where is the cosine of the angle between an arbitrary 

beam and the inward normal to the medium boundary and Tis the optical depth measured 

from the incident flux boundary. The transfer equation is: 

1 2 TT 

‘ ^ pO 1 ,0M f 0) I ,(t /^) cf£d$ + (1^ ) 1 B B 0) 

-I 0 

•where in spherical coordinates P (©) becomes: 

: P(^,0,t,fO=w|l +x [M-ia + (I-i-L 2 ) 1/5 (l-U 2 ) 1 / 2 ** (0-0)]^ (2) 

Substituting (2) into (1) and integrating with respect to 0 yields the transfer 
equation: 

■ u dl ft , ^) = .|( T(|i )H.g.( 1 0 +X .w)l(T, a )dja+.(l-a,)| BB ’- ( 3a ) 

dT z ) ^ 

Evaluating the integral of (3a) via Gaussian quadrature gives: 

U i3? (T ' |J ! ) ‘ S-l(T ' ti i )+ fi-k «>j < ? +XU.H.) I (t ,H.) + (l-a,) l BB (3b) 
where the aj’s are the Gaussian weights, the * s specify the Gaussian directions, and 
k is the order of the Gaussian quadrature. In (3b) a.j = a. and U . = -u ^ . Boundary 
conditions for the transfer equation are given by a specified diffuse intensity field at 
T = 0 and specular reflection at the boundary t = t q given by: 

I(t 0 #H_j) = Rg I(t 0 ,j! .*.)-+ (|-Rg) 1 Bb( T o^ (4) 

where Rg is the reflectance of the medium-substrate interface. 
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For purposes of evaluation of approximate theories discussed below/ right and 

left directed half fluxes and net radiative flux are defined In terms of the transfer 

equation solution field 1( T ) as: 

/I ~ k 

q + (T) = 2rr\ l( T ,U)MdH = 2 tt £ I (t ,U .) 

/ | = 1 I 1 

0 

/-I ~ -1 

q"( T ) = 2rr) l(T,n)ncii = -?tS o { \s.\(r M 
0 i=He M 1 

q (t ) =2rr \ I(t,h)M 1-1 =q + - q~ 


(5a) 


(5b) 

(6) 


Of prime importance in the calculation of the temperature field due to radiative 
absorption is the divergence of the radiative flux which in the plane -parallel case is 
proportional to the derivative dq/dr . In order to avoid numerical differentiation an 
exact expression for dq/dr may be obtained by integration of the transfer equation (3a) 
over the entire solid angle. The resulting expression is: 


■ r - 



It is evident from equations (5), (6) and (7) that the radiative flux and radiative 
flux divergence expressions do not depend explicitly on the anisotropy parameter. Thus 
anisotropy influences those quantities only through the intensity field. Reflectance may be 
determined when the intensity field at the left boundary is known. The hemispherical- 
directional reflectance and hemispherical reflectance are: 

=tt R = <f (°)/q + (°) (8a,b) 

Numerical Solution of the Radiative Transfer Equation 

Replacing the derivative in (3b) by a forward difference quotient forUj>0 and a 
backward difference quotient for M j < 0 yields the set of equations: • 
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f(T+AT^,)-[HAT/i.)]l(T # W i H«AT/2M f ]S a.(l+x*H )\(r,[l ) 

I I 1 

. ; ~ -t(( 1 -tt) ) A T /Uj] 1 BB fo r m > 0 (9a) 

k 

l(T-&T,Uj)=[l-{iT/|Uii JlKT^i^AT^I HiljS^ajd +xH- Uj)l(T,U .) 

• -K(1-U))A T /| Hi I ] I BB for [ij < 0 (9b) 


When an initial approximation fo the radiation field is known and left boundary incident 
. .intensities are specified, equation (9a) can be used fo compute the right directed inten- 
- sities at successive nodal points. When the march from T - 0 to T = T Q is complete, that 
is, when all right -directed intensities have been computed, the rear surface boundary 
--^condition provides values for l( T 0 ,^-j). Equation (9b) is then used to march back to the 
•^incident flux boundary. The initial intensity field approximation is then replaced by this 
■"newly computed field and a second iteration is performed. In practice the iteration 
vrprocedure is continued until two successive iterates for the intensity field differ by less 
-than a specified tolerance. A similar procedure was employed in reference [13] . 

Equations (9a) and 9b) together with the appropriate boundary conditions can be 
written in a matrix vector equation of the form 1 = Mi~*C where I is a vector which has as 
-*&its components the intensities in all directions at all nodes, M is a coefficient matrix, 
r*md C is a known vector determined from the specified boundary data and the values of 
[( l-OJ )A T / 1 J ] I BB . This iterative procedure can be written: 

vj(k+ 1 ) = M|0<) + c for k *= 0,1,2,... (10) 


^rbere !®is the initial approximation to the intensity field. It is well known, [14], that 
%the iterative scheme(lO) is convergent if and only if P (M^T where P (M) denotes the 
^spectral radius of the matrix M. In the case of isotropic scattering If: 


41 J T o ^ 1 ■ 

pi| |^i| (N-l) - 1 - (ajiw/2) 


for all i 


HD 
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•where N is the number of nodes in the finite difference approximation, it can be shown 
that 1 1 Mil < 1 whe re II Mll^ denotes the infinity norm of the matrix M« Since 
p (M)< 1 1 Mil ,[141, if condition (1 1) is satisfied, convergence of the iteration procedure 

(10) is guaranteed. It is interesting to note that condition (11) is satisfied if and only if 
: the coefficient of I(t,^j) in (9) is non-negative. The convergence boundary given by 

(11) is shown in Figure 1 as related to quadrature order. For all cases attempted within 
■this convergence boundary, smooth convergent solutions were obtained. However, it is 
-possible to obtain convergent solutions in cases where condition (1 1) is not satisfied. These 
solutions are characterized by oscillatory Irregularities in intensity near the physical 

-iboundaries and in directions closest to parallel to the boundaries. On the other hand 
■calculations performed fortfJ = .9995 and(&T /min] &jl)greater than about 2.0 diverged 
,*rapidly. For first order Gaussian quadrature, ?t can be shown analytically that the method 
**wfll diverge if (A T /I [X\\ )> (2.0/ w). While only a limited number of runs were made 

for optical thicknesses other than 3.177 and albedos other than .9995, (A T /min |Uj|)> (2.0 /uj) 
appears to be sufficient condition for divergence. It is worthwhile to remark that the 
-^iteration procedure converges more rapidly when newly obtained intensity values are used 
•as soon as they are generated. This Is not suprising since using new intensity values as soon 
*as they are available is analogous to changing from a Jacobi type Iterative procedure to a 
: Gauss-Seidel procedure. •>. 

^In order to verify the accuracy of the iteration technique a computation made with 
-ithe iterative technique was compared with an available unpublished solution obtained by 
--tthe matrix eigenvector technique of [15] . The resulting intensity field using fourth order 
Gaussian quadrature and 21 nodes is shown in Figure 2. The solution corresponds to radiative 
^transfer between two black boundaries through a nonabsorbing Isotropically scattering medium 



(® “ 1.0) where the left boundary and medium are at IIII°K and the right boundary Is at 
555.5° K„ The intensify fields are seen to be in good detailed agreement within about 
1%. This is within the maximum error claimed in reference [15]. For most of the cal- 
culations presented here the order of quadrature was increased to eighth and the number 
of nodes to at least 101. 

Approximate Methods 

The approximate methods for radiative transfer considered stem from two-flux 
models which Have found wide application. Models for isotropic scattering only are 
considered. These models may be related to each other and to the equation of radiative 
transfer discussed above. A link between the transfer equation and the two-flux rela- 


tions known as the Schuster-Schwarzschl Id equations, reference [7], may be established 


by assuming that the intensity field is fully described by two functions of optical depth 
only: 


Kt,u) = 


V(T), 


M > 0 
M < 0 


02 ) 


Here l + represents the value of the directionally independent intensify of rays traversing 
the medium from the incident flux boundary toward the rear ( T increasing), and I" 
similarly represents those rays proceeding forward from the rear (T decreasing). Intro- 
ducing these restrictions and separately integrating the transfer equation over the half 
spaces 0 and \± < 0, the differential equations; 


(l/ 2 )(dl + /dt) =-l + +(w/ 2 )(l + +r) +(!-«) 1 BB (13a) 

-(l/2)(dr/dT) = -I" + (u>/2) (1 + +1") + (1-a) ) l BB (13b) 


are obtained. These are the Schuster-Schwarzschild equations obtained in reference 
[7] in a slightly more general context. The half fluxes and net radiative flux are 
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easily related to the intensities l + and I” using equations (4) and (5): 
q + (T)=2rr ^ I + (t )f-i dU =ttI + (t) 

1 

q“ (t ) = 2n ^ I (t ) |id|i =rr l" (t) 
i' ° 

q - 2tt \ l(T / M>idu =q + (T) -q~( T ) “TT (I + (t) - I “(t )] 

-I 


04a) 

04b) 

04c) 


Using equations (13) and (14) the Schuster-Schwarzschild equations may be expressed 
in terms of the half-fluxes as: 


(dq + /dT ) = ^ -2) q + +iVq~ + 2tt (1-au ) l BB (15a) 

-(dq‘/dr) = fc)-2) q“ +uj q" f " + 2rr (l-a>) l &B (15b) 

Equations (15) are the isotropic form of the transport equations given by Bergquam and 
Seban [11]. In reference [1 1], equations (15) together with an energy equation were 
solved iteratively on a digital computer to evaluate the two flux method for emission- 
coupled cases. 


The albedo ^ and optical thickness r have similar interpretations in the equation 
of radiative transfer and in the Schuster-Schwarzschild equations. It is well known that 
these quantitifes may be related to constant scattering, absorption and extinction co- 
efficients S, K and p through 

S 

T - py, p - S + K, u) - s + £ (16a, b,c) 

where y Is the physical depth measured from the incident flux boundary. It is through 
these coefficients that direct associations may be made with another two-flux model to 
be discussed later. Transforming the independent variable in equations (15) to y and 
introducing the scattering and absorption coefficients ? the Schuster-Schwarzschild equations 


become : 
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(dq+/dy) = -(2K+S)q + + Sq“+2n KI bb (17q) 

- (dq-/dy) = -(2K + S ) q- + S q + -t^T K l BB (17b) 

An entirely separate tradition in two-flux radiation field analysis has developed 
im the paint and paper industry,, The basic equations attributed to Kubelka and Munk 
.[8] may also be formulated in terms of half fluxes q + and q ~ and the physical depth y. 
rThese equations, however, include independently defined scattering and absorption 
^coefficients s and k. The equations are formulated by requiring that each half flux be 
augmented by scattering from the opposing half flux and by emission, and be diminished 
by scattering to the opposing half flux and by absorption. The Kubelka-Munk 
-differential equations are: 

(dq + /dy)= ~(s + k) q + + sq +tt k , (18a) 

(dq“/dy) = (s +k) q“ - sq + -TT k IbB (18b) 

* Comparison of equations (17) and (18) shows that the simple relations: 

s = S, k =2K (19a, b) 

'allow association of the Kubelka-Munk equations with the Schuster-Schwarzschild 

equations and hence with the equation of radiative transport. The optica! depth and 

ualbedo may then be related to the Kubelka-Munk coefficients through: 

2s u 

W = 2TTk' T=(s+ 2 )> ' (20a,b) 

/ ’ • 

"While the association of the Kubelka-Munk equations with the Schuster- 
^Schwarzschild and radiative transfer equations discussed above is satisfying in its 
simplicity, it is clear that the relations (19) and 20) are at best approximations 
■ywHich will be most realistic when the intensity distribution in the medium approaches 
the distribution given by equation (12). Many alternatives to the equations (19) are 
^possible. Klier [10] has, for instance, shown that the Kubelka-Munk equations are 
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formally identical to a class of solutions of the radiative transfer equation for isotropically, 
highly scattering media. Under this approach the Kubelka-Munk s and k may be related 
to the transfer equation parameters through 

s=SF ] (|), k = Kf 2 (|) (21 a,b) 

The nonlinear functions f^ and f 2 are tabulated in [10]- For weak absorption the reference 
gives f j - .75 and ^ ~ 2-0. Thus, in this case, equation (19) differs from the Klier 
equations only in the constant coefficient of the scattering coefficient relation. 

In the Results and Discussion section, Kubelka-Munk solutions for reflectance, 
radiative flux and radiative flux divergence are compared with transfer equation solutions 
using both relations (19) and (21). These solutions are given below: 

Exact Kubelka-Munk Solution for Finite Thickness 

Following Hamaker [8] the exact solution of the Kubelka-Munk equations (with 
negligible emission) for the half fluxes may be written: 

q + =A(l^)e° y + B. (1 +Y )."% . q"=*A(l +-Y )e^ + B (1-Y ^ (22 a,b) 


where 


a =7 k(k+2s) , Y =7 k/(k+2s) = a /(k+2s) 


and A and B are determined by the boundary conditions. 


(23 a, b) 


* .For the c ase o f a scatteri ng , abso rbi ng, non-e mitti ng mediumwi th inte nse , 
diffuse incident flux q + (0) and opaque non-emitting substrate with reflectance Rg, 


the coefficients A and B are: 

. _ q + (0)e^ 5 [y(l+RBH1-RB)] 

2 £[y 2(1 +Rb)"K 1 “Rg )] sin h (a 6 ) _ + 2 Y cos h (ffS )} 

B ^ q + (0)e gS [Y(l-m R )-Kl-R R )] 

2 {[y^ 0 — Rg)l sinh (06 ) + 2Y cosh p 6 )} 


(24a) 

(24b) 


where 6 is the thickness of the scattering medium. 


(27) 

(28) 


This solution is used together with equations (14) to provide an initial radiation field for 
the numerical iteration of the transfer equation discussed earlier. The Kubelka-Munk net 
radiative flux and radiative flux divergence may be found by using equation (14c) and 
eliminating derivatives using the Kubelka-Munk differential equations (18): 

W-T 1251 

■ £-£ + -£'--^o 

The hemispherical reflectance of a layer of infinite thickness, , and of a 
finite layer may be obtained from (8) and (22) as: 

p <CM. = III = 1 j —(—+2) 

q^W s J s S * 

and 

(1/IU(Rr-IU - m (Rn-lA.) exp [s6 (l/Rco-R,)] 

R “ (Rg-Ros) - (Rb“VR®) ex P t s6 0/R»~R®)] 

Approximate Solution for Weakly Absorbing Media 

A good volume reflector is characterized by low absorption coefficient and high 
scattering coefficient. In this situation an existing approximate solution of the Kubelka- 
Mjnk equations with negligible emission may be useful. This approach, used in [9], 
involves solution of the Kubelka-Munk equations for the case of zero absorption coeffi- 
cient subject to a specified incident diffuse flux and specified overall reflectance . The 
resulting solutions may be written as: 

q*(y) = q + (°) n-syO-R)l/ q“(y) = q + (°) [R-syO-RM < 29 a ' b ) 

This approximation implies that the net radiative flux is given by: 

q = q+- q -=q + (0).[l-R] < 30 > . 

It is evident that in this approximation q is independent of depth. Using equation (26) 
the radiative flux divergence may be approximated by combination with equation (29): 

(dq/dy) — -k q + (0) [1 + R - 2 sy ( 1 —R)] < 31 > 
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- In these approximations the reflectance R is evaluated using equation (28). Other suitable 
approximations may be used however. Note that expressions (29) and (31) are clearly limited 
on physical grounds to depths for which the terms in brackets are positive. 

RESULTS AND DISCUSSION 

When an intense, diffuse radiative flux is incident on a highly scattering and weakly 
absorbing medium it is well known that an intense, diffuse reflection occurs. The reflected 
flux, for medium indices of refraction close to one, originates as a result of penetration in 
depth of the incident rays, rather than as a surface phenomenon. Internally, multiple scat- 
tering processes cause a redirection of a large fraction of the incident radiant energy back 
toward its surface of entry. The ultimate re-emergence of the majority of the incident and 
scattered photons is made possible by the absence of strong absorption within the medium. 

The manner in which energy is redistributed is shown In Figure 2. There the Intensity 
-variation with depth Is presented for a nonabsorbing, isotropically scattering medium, (yj=1.0), 
with a black substrate. It is evident that the Incident transmitted rays diminish in intensity 
with depth even in the absence of radiative absorption. This effect Is more pronounced as 
the path length increases ii- approaches zero). Photons penetrating to the black substrate 
are absorbed there and new rays emerge directed toward the incident flux boundary. These 

rays, initially equal in Intensity, are augmented in intensify as they progress towards = 0 

/ 

*’ 4 

due to scattering from each of the other rays. Again rays traversing greater path lengths 
ore influenced more by scattering than those closer to the surface normal. 

The influence of the substrate reflectance on the dimensionless intensity field is 
shown in Figure 3 fora weakly absorbing, isotropically scattering medium with negligible 
emission. The lines and symbols show Intensity distributions for substrate specular reflec- 
tances of 0 and 0.8 respectively. Comparison of the two intensity families shows clearly 
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the importance of substrate reflectance even for optical thicknesses of the order of three. 

The distributions of dimensionless intensity at the Incident flux boundary (*1 = 0) are also 
hemispherical-directional reflectance distributions. The monotonic decrease in reflectance 
toward the outward normal to the incident flux surface Is apparent. Evidently the increased 
scattering influence due to path length discussed above is controlling this distribution. 

The influence of quadrature order on the intensity field is demonstrated in Figure 3 
also. It is evident that quadrature orders above four are necessary only for filling in the 
distribution, as has been pointed out elsewhere in a related problem [4]. 

It is indicated above that the Influence of emission from the medium and substrate 
may be negligible under an intense incident flux. For the incident field considered in 
Figure 3, corresponding to a black emitter at 3000°K, calculations indicate that the 
intensity field changes by less than 1.5% when the medium and substrate temperatures 
are changed from 30Q°K to 1000°K. Thus neglect of emission is well justified for a large 
class of dielectric solids at moderate temperatures. 

The influence of anisotropic scattering on the internal intensity field is presented 
in Figure 4 by comparing transfer equation solutions for strong forward scattering (x = 1.0) 
and strong backward scattering (x = -1.0) with the Isotropic solution (Rg = .8) of Figure 3. 

It is clear that even with this extreme anisotropy the basic character of the distributions is 
unchanged. It may be observed in Figure 4 that forward scattering reduces the attenuation 
of incoming rays and that this effect becomes more pronounced with increasing depth. 

While the substrate -re flee fed Intensity for forward scattering is greater than for the isotropic 
case, reduced augmentation by scattering into rays directed toward the incident flux 
boundary produces a lower Intensity level emerging from the medium than for Isotropic 
scattering. Thus the medium reflectance is reduced due to forward scattering with respect 
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to the isotropic case. This may be viewed simplistically as increased trapping of energy 
due to deeper penetration by forward scattering and due to absorption of this energy in 
the medium and at the substrate. Conversely for backward scattering less energy pene- 
trates to a given depth to be trapped than for isotropic scattering . As a consequence, more 
energy is returned to the incident flux boundary in the emerging rays, leading to Increased 
reflectance for backward scattering media. Hemispherical reflectances calculated for 
each of these cases are tabulated on the figure. 

The total diffuse reflectance of isotropically scattering volume reflectors is con- 
sidered in Figure 5 as a function of albedo, optical thickness and substrate reflectance. 

Values computed from the Kubelka-Munk reflectance equation (28) show the transition 
from optically thin to optically thick media for several values of albedo and substrate 
reflectance. The solid tines show the Kubelka-Munk reflectance interpreted through the 

V 

relations (19) while the broken lines utilize the Klier relations (21) and the tabular values 
of f] and f 2 from [10]. It is evident that the difference between the two representations 
grows with decreasing albedo. 

Values of hemispherical reflectance based on solutions of the equation of radiative 

transfer are also presented in Figure 5. In addition to values generated in the present study, 

published results from Ui and Ozisik [16] for media of finite optical thickness and Giovanelli 

[17] for semi-infinite media are included. Three of the Ui and Ozisik cases were run for 

- comparison purposes with the present iterative program. Two of these were in agreement 

to three significant figures using 8th order quadrature and 101 nodes. The third, within 1.4% 

of the LI? and Ozisik value, (r = 5.0 and tf) = .95) was rerun for 201 nodes and the results 

o 

were extrapolated in terms of the reciprocal of the number of nodes. Extrapolation in this 

i 

way to an infinite number of nodes yields agreement with the Ui and Ozisik value in the 


fourth significant figure. 
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The transfer equation solutions must of course be regarded as correct values for the 
reflectance of the scattering, absorbing medium. The two representations of the Kubelka- 
Munk solutions both show the same qualitative trends as the transfer equation solutions. 
Quantitatively, the Kubelka-Munk equations using the relations (21) show better agreement 
with transfer equation solutions than when equation (19) is used. 

Flux and Flux Divergence 

The radiative flux and radiative flux divergence variations with depth are presented 
in Figures 6 and 7 for the intensity field of Figure 3 (R B = .8). In these figures both repre- 
sentations of the Kubelka-Munk exact and approximate fluxes and divergence are compared 
with solutions of the transfer equation. The Kubelka-Munk results are seen to be reasonable 
approximations particularly when the Klier representation [10] is employed. For this partic- 
ular case (T 0 =3.177) the Kubelka-Munk exact and approximate values are equally good. 
This is fortuitous except near the incident flux boundary. It can be shown that approximate 
and exact Kubelka-Munk expressions are identical at r\ = 0 and tend to deviate increasingly 
with depth. In all calculations performed, including several not presented here, Kubelka- 
Munk divergence values were in excellent agreement with transfer equation values at the 
incident flux boundary. In view of the limitation of the approximate solutions discussed 
following equation (31) two additional optical thicknesses are shown in Figure 7. It is evi- 
dent that both types of Kubelka-Munk solution deteriorate at the highest optical thickness — 
.the approximate solution the more rapidly. 

. Also shown in Figures 6 and 7 are flux and flux-divergence distributions for the 
anisotropic cases of Figure 4. The Increased penetration of the incident flux and net flux 
for forward scattering and decreased penetration for backward scattering as discussed above 
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are evident in Figure 6. Increased energy loss from the radiation field for forward 
scattering and decreased loss for backward scattering are evident in the corresponding 
increase and decrease in the divergence distributions in depth in Figure 7* 

The influence of the albedo of highly scattering layers on the net radiative flux 
and on a flux divergence parameter is shown in Figures 8 and 9. It is evident that the 
Kubelka-Munk exact solutions provide good approximations over the entire range of w 
considered, particularly when equation (21) is employed. It is also apparent that the 
Kubelka-Munk approximations are useful only for > .995. That this is not a fatal 
limitation is evident from Figure 5 which shows that this range Includes optically thick 
volume reflectors with reflectances as low as ,85. The limitation of the approximate 
divergence equation (31) discussed earlier is apparent in Figure 9. Notwithstanding 
their limitations, the Kubelka-Munk approximations may be utilized in the study of high 
performance volume reflecting heat shields. 

CONCLUDING REMARKS 

Accurate numerical solutions of the equation of radiative transfer have been 
developed for the evaluation of approximation methods for reflectance and radiative flux 
parameters. Convergence and divergence criteria are given for the numerical solution of 
the system of coupled differential equations used. New hemispherical reflectance values 
are presented and compared with existing values In the literature. The Kubelka-Munk and 
Schuster Schwarzschild two -flux theories are obtained as consequences of the equation of 
radiative transfer under the assumption of directionally independent intensity fields. This 
approach yields simple linear relations between transfer equation and Kubelka-Munk 
scattering and absorption coefficients. An alternate interpretation of the Kubelka-Munk 
- theory as an exact theory of an asymptotic solution of the equation of radiative transfer 
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yields slightly different nonlinear relations between the coefficients. Exact and approximate 
Kubelka-Munk solutions were compared with exact solutions of the equation of radiative 
transfer using both interpretations. The exact Kubelka-Munk solutions are found to provide 
realistic approximations to transfer equation results over a wide range of conditions while 
approximate Kubelka-Munk relations are applicable over a more restricted range of conditions. 
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Figure 5: The Hemispherical Reflectance of Isotropic Scattering Layers Mounted on a 
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Figure 6: Comparison of Radiative Flux Distributions Calculated from the Transfer 

Equation and from the Kubelka-Munk Theory. T =3.1 77, uj = .9995, 

Rg = .8, Negligible Emission. 

Figure 7: Comparison of Radiative Flux Divergence Distributions Calculated from 

the Transfer Equation and the Kubelka-Munk Theory. 05= .9995, Rg = .8, 
Negligible Emission. 

Figure 8: The Influence of Albedo on the Net Radiative Flux at Three Optical Depth 

T 0 = 3. 177, Rq = .8, x = 0, Negligible Emission. 

Figure 9: The Influence of Albedo on the Radiative Flux Divergence at Three Optica 

Depths. =3.177, Rg = .8, x = 0, Negligible Emission. 
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